INTERNAL SPACE RENORMALIZATION GROUP METHODS FOR ATOMIC 

AND CONDENSED MATTER PHYSICS 



Janos Polonyi 

Laboratory of Theoretical Physics, Louis Pasteur University, Strasbourg, France, 
Department of Atomic Physics, Lorand Eotvos University, Budapest, Hungary 

The functional renormalization group method is used to take into account the vacuum polar- 
ization around localized bound states generated by external potential. The application to Atomic 
Physics leads to improved Hartree-Fock and Kohn-Sham equations in a systematic manner within 
the framework of the Density Functional Theory. Another application to Condensed Matter Physics 
consists of an algorithm to compute quenched averages with or without Coulomb interaction in a 
non-perturbative manner. 



I. INTRODUCTION 

The renormalization group will be used in this talkQ as an algorithm to solve strongly coupled quantum field theories 
without any intention to gain insight into the scale dependence of the dynamics. We shall study the system of non- 
relativistic electrons propagating in the presence of static external potential. We consider non-relativistic systems 
because the comparison with experiment is more direct than in the relativistic domain. The common challenge in 
both regimes is to trace the polarization effects in the vacuum, i.e. the particle-hole or the particle-anti particle 
fluctuations in the non-relativistic or relativistic region, respectively. 

We shall consider two different cases. First, the electrons will be placed in an external localized field and one finds 
a problem with inhomogeneous ground state, characteristic of Atomic Physics . Second, we assume the presence of a 
static random external impurity potential which obeys a Gaussian probability distribution and we introduce a scheme 
to compute the quenched averages of kinetic transport coefficients in a translation invariant manner. The common 
technics applied is a generalization of the functional renormalization group |l[ where the role of the running cut-off is 
played by an arbitrary control parameter which generates differentiable changes in the dynamics. 



II. RENORMALIZATION IN THE INTERNAL SPACE 



The degrees of freedom are eliminated successively in the renormalization group method, for instance the Kadanoff- 
Wilson blocking strategy orders the modes according to their scales in the external space, i.e. in the the space-time 
(energy-momentum). One may construct a blocking procedure where the blocking proceeds in an order determined 
by a scale in the internal space, the space of the field amplitude. A well known example is the Callan-Symanzik 
equation where the fluctuations with larger and larger amplitudes arc taken into account as the mass of the particles 
is lowered. 

It is easy to generalize this scheme for any parameter in the dynamics by means of functional techniques Qj. In 
fact, let us suppose that the dynamics contains a control parameter A and the generator functional for the connected 
Green functions is given by 



where / • g = J dxf x g x . The evolution equation 
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is valid without assuming the existence of any small parameter. The way this equation is obtanied is reminescent of 
the derivation of the Schwinger-Dyson equations except that the latter contains the complete action and the former 
the suppression controling part, d\S\, only. It is worthwhile noting that this relation appears to be such a functional 
generalization of the Hellman-Feynman theorem || , 

d\E\ = (E x \d x H x \E x ), E X \E X ) = H X \E X ), (3) 

for time dependent processes and general matrix elements which remains compatible with approximations one employs 
to solve the functional differential equation (^) in restricted functional spaces. 

In the examples presented here the electric charge and the average strength of the impurity field fluctuations will 
be chosen as control parameters. The corresponding evolution equations will be written for a local functional, the 
effective action for the density and the current. The effective action of the physical system can be obtained by 
integrating these equations from an artificial perturbative initial condition imposed at weak Coulomb interaction and 
disorder into the physical regime. 

Similar schemes could be constructed by means of the more traditional renormalization group method, based on 
the external space. But the present version has the following two advantages. First, it preserves gauge invariance, 
an important feature in the computation of the electric conductivity. Second, it avoids the artificial discontinuities 
during the evolution which may occur when the saddle point structure is changed by the running cut-off . 



III. DENSITY FUNCTIONAL THEORY 



The density functional theory [|5|J6|] is a powerful method to describe bound states induced by external potential. 
The quantity of central importance is the density functional, -E gr [p], the ground state energy when the the density 
is constrained to be p. According to the the Hohenberg-Kohn theorems the ground state energy of the original 
problem is the minimum of the density functional, E gI = E gI [p gr ] . The shortcoming of the traditional approach is 
its phenomenological nature, the lack of a constructive definition of the density functional. We propose that the 
effective action T[p] for the density p 0] provides a clear and generally applicable definition of the density functional. 
Furthermore, the internal space renormalization group method leads to a systematical approximation scheme which 
turns the phenomenological knowledge about correlations into an improvement of the approximation by choosing the 
ansatz for T[p] in an appropriate manner 

The generator functional for the connected Green function is defined as 



where j a ,x is the electric density (a — 0) and current (a — 1,2,3). We shall consider the density Green functions 
and use a = only in this section. The Legendre transform of the IFfcr], the effective action T[p] for the density p 
reaches its minimum at the ground state density, f3E gI = r[ ( o gr ], and will serve as the density functional. The effective 
action is first obtained for imaginary time, at finite temperature, T = 1//3, and the projection onto the ground state 
is achieved by taking the zero temperature limit f3 — > oo. After the rescaling e — > Ae of the electric charge one can 
derive the evolution equation 



d x T x [p] = -Xe 2 | p-^-p + Tr 



SpSp 
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Wyi,t^,t%^y,t)d x —, (5) 

x,y,t »7T|x-y| 

where tp denotes the electron field operator and J = J dx. The second equation shows that the evolution equation 
simply follows the change of the Coulomb energy according to the Hellman-Feynman theorem. 

The initial condition for the effective action is imposed at weak Coulomb interaction, Ao ~ 0, where we find in the 
leading order perturbation expansion for spinless electrons 

m = \{p -P*)-o- l -(p-P*)-P~-P +c[ P ] (6) 

with C[p] = -TrlogG- 1 + \Tr\ogD- 1 + 0{e 2 ) + 0{p A ). Here 
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= S*,« (d t - £-A x - M + V^) (7) 

is the propagator for non-interacting electrons in the presence of an external potential V and chemical potential 
Mi G X . XJ = —G x ^yG VtX stands for the particle-hole propagator, Q = [G _1 + e 2 /A]~ 1 is its improved one-lop version, 
D = [—A + e 2 G}~ 1 denotes the photon propagator and p* = —G x _ x . By minimizing the effective action one recovers 
the Hartree-Fock energy functional in 0(e 2 ). The non-classical nature of the exchange contribution to the interaction 
is reflected in the fact that it arises from the photon fluctuation determinant, the term TrlogZ? -1 . 

The evolution equation (^|) should be projected into a restricted functional space in order to make it more man- 
ageable. For this end we introduce the multi-local truncation scheme for the effective action. The free fc-local cluster 
approximation, corresponds to the functional form 

r M = E II E / <'-o IT • ■ • pZ (8) 

3=0 j=l n 3 J 

with arbitrary The effective action in the constrained fc-local cluster approximation, Cfc, can be written 

in the same manner except that the functions r™* are parameterized. The interactive electron field operator is 
assumed to be of the form 

oo 

1pt,x = E c n,t^n,m, (9) 
n=l 

and the propagator will be written as 

N oo 

G x , x , = J2 e ~ EAt ~ t ' )lS '^n, x '+ E S«.n'(*-*')*n,x*;' > x' > (10) 
n— 1 n,n' — 1 

for t < t' leaving ^ n:X , E n and g n ,n'(t) as parameters. 

We introduce a local density-dependent self-energy, a(x,p x ), in the photon propagator and assume the form 



C[p] = -Tr log G- 1 + ^Trlog^- 1 + a] + j U(x, p x ) + £ p n ■ 7 <" 



rn ) m 



p m (ii) 



which is a ji ansatz with infinitely many higher order constrained clusters arising from the photon fluctuation 
determinant. The corresponding evolution equation, considered at the minimum of the effective action, 

can be written in a variational form as 

6H[V*,V,E,p] 6H[V*,y,E,p] dH[^*,^,E,p] 



5 *n,x \ P=Psr \ p = Psr 9E n lp=Psi 



0, (13) 



after the appropriate choice of a constant Tq. A generalized HF functional 

*,#,£,/>] = H iT [^*,^,E} +H c [^*,y] +H ph [**,*,p] (14) 
was introduced here which is the sum of the one-particle, exchange, direct and interaction pieces, 

" r r / *2 



^ h [**, p] = A Tr log (-A + a(p) + X 2 e 2 d) , 

\2 2 W /• 1 
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The variational form (13) has the following remarkable features. The generalized Hartree-Fock functional includes 
higher order radiative corrections which involve time-dependent, dynamical quantities. Furthermore, it shows that 
the evolution, imposed at the ground state where the effective action is the best approximated, automatically involves 
an optimization with respect to the choice of the quasi-particles, the wave functions ^ n in the electron field operator 

®- • 

We shall consider two simple approximation schemes for the evolution equation (|13|): In a C2 truncation we keep 
the 0(e 2 ) and 0(p 2 ) perturbative effective action and recover the usual Hartree-Fock equations for the single particle 
wave functions and energies E n . Another f\Ci approximation scheme is where one keeps the local potential U(p) 
arbitrary but sets 7 = 0. It is important to realize that the transformation 

2 

U{x,p) -> U(x,p) +T)(x,p), <r(x,p) -> a(x,p) + — — i){x,p) (16) 

leaves the effective action invariant. As pointed out above the photon fluctuation determinant gives the exchange 
contribution and is density independent. But the transformation ( |l6| ) can be used to trade a possible, non-perturbative 
density-dependent photon self-energy term a into a local potential U. The result is a generalization of the Kohn-Sham 
scheme pj, namely Eqs. (|l^) applied for the functional 

n KS [*• , *, e, p] = n b [**,*, e] + n° s [** , *, P ] + ?4 S p], (i7) 

where 

N 

n^*^,p}^-n s X 2 e 2 E f ~Ax 4 , , Py,o (18) 
1 A 2 e 2 A 2 e 2 

K KS [^, tt, p] = --p ■ ■ G ■ ■ p. (19) 

The ground state energy is determined by the differential equation 

r -11 A 2 e 2 - 2Ae 2 - A 2 e 2 

d, J U(x, 0) = -Ae 2 Tr G-^ - ^p* ■ ■ G ■ ■ G ■ ■ p* (20) 

together with the initial condition U — and contains the exchange term. 

IV. QUENCHED AVERAGES 

The impurities represent a challenge in Condensed Matter Physics since they appear to be static from the point 
of view of the measurements and are distributed randomly. The usual way to take them into account is to average 
the connected Green functions, the logarithm of the partition function over the impurity distributions j9j. This 
average can be obtained either by analytical continuation in the number of replicas [lFJ, or by the introduction of 
fictious particles related to the real ones by super-transformations ]TT[ | or by using the Keldysh contour in computing 
loop-integrals fl^] . The functional renormalization group idea offers an alternative algorithm to compute quenched 
averages which goes further than these methods in being fully non-perturbative and in allowing annealed interactions 
fl3f . For the sake of simplicity we constrain the present discussion to the non- interacting case only. 

Consider the generator functional for the Green functions of the density and current, introduced in Eq. (Q), and 
write its quenched average as 

w[a] = J^ J '^;+^^1, (21) 

The Legendre transform of W[a], the effective action T[p], satisfies the evolution equation 



dxT[p] 



- 9 -l 
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The gradient expansion ansatz for the functionals T[p] seems natural when the kinetic transport coefficient are 
sought, 




d u T tt id - -dlT u dl - 5 g2 r«A 



Po,. 



^o,,(r ts + d„r ts id )d kPk , z - ^ M (r st + d^r st ido)d kPQ ,. 



1 



Sj, k ( r ss + d u r ss id - lo 2 r ss d 2 



-r T («F' fc A - djd k ) - ria,-s fc 



p k , z + y tss ~ Ufa) 



(23) 



The r coefhcients are functions of the control parameter A and the density po and have been computed in the leading 
order of the perturbation expansion in g. Such an approximation is reliable for weak disorder, k = gm /2Trh 3 p F < 1, 
and one finds the form P(n)/p 7 p for the different T-functions where P(k) is simple regular function, n > and r — > 
when k — ► oo. Some combinations of the T-functions appear in the Kubo formula, e.g. the electric conductivity, 
<?j,k = ^,fee 2 a w r 5 V(r ss ) 2 , the diffusion constant, D = d q *T u / \T U ) 2 , and the quantities T ss and T u = -d 2 pQ U give 
the density and current susceptibilities, respectively 

When the vertex corrections are ignored then the evolution equation can be written as 

a A r = o for r = a w r",^r ss , 

d x r = gd 2 P0 rJ^ for r = r ss ,r T ,r L , 

dxl = g^d P0 + 1) ( d 2 pnl f ±) for 7 = 9 2 7 ", d q * 7 tt 



9a7 = .9 



f 



f 



P 2d 2 y 



(24) 



where a 2 7 " = (po^ + l)5 2 r«, d q2 j tt = (p a d po + l)d q2 T tt , 7 * s = {p Q d pa + 2)r ts , 9 2 7 ts - (p 5 P0 + 2)9 2 r ts , and 



d P = P 2 ( Po d po + i)d q2 r u + r" + - [(p 5 po + i)r ts - r st ] " 



P 2 r L + r s 



(25) 



It is easy to see that the integration of these equations resums all correlation insertion of the particle-hole loop diagram 
in the given truncation of the gradient expansion, in particular the maximally crossed diagrams which are responsible 
for the cooperon pole needed for non- vanishing conductivity p"4| ] . 

The terms O(d 2 ) of the evolution equations correspond to a formal diffusion process on the plane (po, <?) considered 
as space-time. One can verify that the remaining terms which contain pad po generate similar spread, as well. The 
diffusion constant in the po _s pace, the integral on the right hand sides, is small at the initial condition for strong 
disorder and becomes large for weak disorder. Therefore the T-functions will approach zero for weak disorder due to 
the strong diffusion. This is supposed to generate weak localization. 

The localization transition should occur when some or all 0(d 2 ) terms of the effective action are vanishing. In fact, 
according to the Kubo formulae the mobility is proportional to certain 0(d 2 ) terms in the effective action. Another 
reasoning is to recall the reduction formulae for a single particle which gives the scattering amplitudes in terms of the 
residuum of the connected propagator on the mass shell. It remains to be seen by detailed numerical studies of the 
evolution equation in 3 dimensions whether the po-dependence of the integrals on the right hand sides provides such 
a self-acceleration of the diffusion process in the po _ space which cancels some T-functions at finite value of g. 

We note finally that there is a formal analogy between localization and the phenomenon of the quark confinement. 
According to the haaron-model of the QCD vacuum the quark propagator is 0(p~ 4 ) and the quark confinement 
appears as a localization in the space-time |15|] . 
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V. SUMMARY 



Two applications of the internal space renormalization group were given demonstrating the possibility of a new, 
systematical non-perturbative method to tackle the bound state problem at least when the bound state formation is 
triggered by an external potential. 

This method is in its infancy and the outline of its formal structure is sketched only. But we believe that it 
can provide an accuracy and flexibility superior to other procedures when formal computer algebra and numerical 
integration are combined in deriving and solving the evolution equation in a sufficiently rich functional space. 
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